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Abstract. We show that beta ensembles in Random Matrix Theory with generic real 
analytic potential have the asymptotic equipartition property. In addition, we prove a 
Central Limit Theorem for the density of the eigenvalues of these ensembles. 



1. Introduction 

In this paper we study asymptotic properties of the density functions of certain mea- 
sures known as beta ensembles that arise in random matrix theory. Namely, we consider 
probability distributions in M. N of the form 



(1) *#(Ai, A 2 , . . . , X N ) = J] |A, - A/ n 



N 

l3NV(\i)/2 



Zn(P) . 



e 



where the potential V is a real valued Holder function satisfying the growth condition 
(2) 7(A) > 2(1 + e)log(l + |A|) 

for all sufficiently large A. 

Beta ensembles have attracted significant interest in recent years (see e.g. [8], [9], [12] . 
P3], P3], [IS], [21], [25], [27], [2H], [29], [30], [32], and references therein). Below, we briefly 
mention four classical beta ensembles, namely Hermite (Gaussian) beta ensemble, Circular 
beta ensemble, Laguerre (Wishart) beta ensemble, and Jacobi (MANOVA) beta ensemble. 

• Hermite beta ensemble: Let Ai, . . . , Xn £ M be random variables with joint density 
function with respect to the Lebesgue measure given by 

1 N 

0) ps^k a 2 , • • • , a.) = n - ri e "^- 

For (3 = 1,2, and 4, the distribution (EJ is known as the joint distribution of the eigen- 
values of a random matrix from the Gaussian Orthogonal Ensemble (GOE), Gaussian 
Unitary Ensemble (GUE), and Gaussian Symplectic Ensemble (GSE) respectively. 
Dumitriu and Edelman (see [12], [3]) introduced tridiagonal real symmetric random 
matrices with i.i.d. centered Gaussian random variables on the diagonal and x dis- 
tributed independent random variables on the upper and lower sub-diagonals such 
that the joint distribution of the eigenvalues is given by (J3]) for arbitrary (3 > 0. 

The next three classical ensembles do not formally belong to the class ([I]) since the 
particles are distributed, respectively, on the unit circle, positive half-line, and the 
interval [—1,1]. 

l 
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• Circular beta ensemble: Let Ai,...,Ajv £ [0, 2tt] be random variables with joint 
density function with respect to the Lebesgue measure given by 

(4) P c /(Xi, A 2 , . . . , A,v) = -i- J] \e^ - e<* 

Z n(P) k<j 

• Laguerre ensemble: Let Ai, . . . , Xn £ [0, oo) be random variables with joint density 
function with respect to the Lebesgue measure given by 

1 N 

(5) ^(Ai, A 2 , . . . , A.) = -±- J] [A, - A/ J] Ar X e-^, 

Z JV \P) i<j j=l 

where a > 0. 

• Jacobi ensemble: Let Ai, . . . , A^ G [—1,1] be random variables with joint density 
function with respect to the Lebesgue measure given by 

1 N 

(6) p^(a 15 a 2 , . . . , a*) = -r n \\ - Aii" na - A.r-^i + A,r \ 

where fi, v > 0. 

As in the Hermite case, the values (3 = 1,2, and 4 (with a = |(n — m + 1) — 1 in the 
Laguerre case and /i = ^(n\ — m + 1) — 1, /x = f (ra 2 — m + 1) — 1 in the Jacobi case), 
correspond to classical ensembles of random matrices (see e.g. [20], [12]). In the Laguerre 
case, the random matrix ensemble for arbitrary positive /3 was constructed in [12]. For Jacobi 
and Circular beta ensembles, three- and five-diagonal matrix models were derived by Killip 
and Nenciu in |19j . 

It is well known (see [7], [18], [11]) that if the potential V in (pQ) is sufficiently smooth 
(e.g. when V is Holder continuous), then there exists the equilibrium measure fi v which is 
absolutely continuous with respect to the Lebesgue measure and has compact support. Let 
us denote its density by p v . The equilibrium measure maximizes the functional 

(7) S v (ji) := J J log \x - y\dfi(x)dfi(y) - J V(x)dfi(x) 

over the space of the probability measure on M. Note that £y(fi), up to a factor — 1, coincides 
with the energy functional. We denote by £[V] := £y(fj, v ) the maximum value of the 
functional attained at the equilibrium measure fi v . The marginal density of the ensemble 

Pi? (A) = / p l3 N {X,X 2 ,...,XN)dX 2 ...dX N , 

JRN-I 

weakly converges to the equilibrium density p v . In addition, the support of fi v and the 
density p v are uniquely determined by the Euler-Lagrange variational equations 



(8) 2 J \og\x -y\dp v {y) - V{x) < /, iGl, 

(9) » = p"(»)>a 

In the Hermite case (V(x) = x 2 /2), the equilibrium measure is the Wigner semicircular 
distribution p sc (dx) = ^V4 — x 2 l[^ 2 ,2](x)dx. 
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In addition to t^j, we assume throughout this paper that the potential V is real analytic 
and generic (see Theorem 12.11 for details). In particular, the support of the equilibrium 
measure consists of finitely many intervals (see e.g. [2j). We prove that the logarithm of 
the density, after appropriate normalization, converges to a constant almost surely. This 
is called the asymptotic equipartition property. The precise statement is given in Theorem 
12.11 By analogy with the Shannon-McMillan-Breiman theorem, the limiting constant can 
be interpreted as the entropy of the corresponding measure. A corollary of this result states 
that in the limit when the size of the considered matrices goes to infinity, there is a set of 
measure almost 1 such that at all the points of this set the density is almost the same. 

Analogous results for the Plancherel measure were conjectured to be true by Vershik and 
Kerov [31] and proved in [10] . For a one parameter deformation of the Plancherel measure, 
called the Schur-Weyl measures, equivalent results were obtained in [211 122] ■ Note, that in 
the case of the Plancherel and Schur-Weyl measures it is only known that the normalized 
logarithm of the density converges to a constant in probability. The question of almost sure 
convergence for those measures is still open. 

In addition to the asymptotic equipartition property we prove Central Limit Theorem 
for the logarithm of the density (see Theorem I2.3p . For general random matrix models 
Borodin and Serfaty have studied asymptotic properties of a similar statistic, which they 
call "renormalized energy" [6]. They explicitly calculate the limit of the expectation of the 
"renormalized energy" in the case of /3-sine processes for /3 = 1, 2, 4 and some 2-dimensional 
point processes. The result of Theorem 12.31 in the special case of the circular beta ensemble 
appears in [5]. Central Limit Theorem for random variables of the form ^^(Aj), where 
Aj are the eigenvalues and g is a bounded continuous function, was obtained by Johansson 
for the circular ensemble [T7] and later for ensembles with more general potentials JTS] (see 
also [23J for results on non-Gaussian limiting fluctuation of linear statistics in the multi-cut 
case). A large deviation principle for spectral measures of certain classes of beta ensembles 
of random matrices was established in [5] . 

For the classical beta ensembles (Hermite, Circular, Laguerre, and Jacobi beta ensembles) 
our results follow from the Selberg Integral [26] since it allows us to obtain explicit formulas 
for the Laplace transform of the normalized logarithm of our density. We obtain the pointwise 
convergence and the Central Limit Theorem from the asymptotic of the exponential moment. 
In the general case, we use recent results [27], [7J, [2B] about the asymptotic expansion of 
the partition function of ([Q). 
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2. Main results 

2.1. Asymptotic equipartition property. Consider a probability measure on M. N with 
density function P^ defined in (CQ). Consider the (infinite) product of these probability spaces 



p JV 

uuuaumy mcaouic un a& . — ir\ a m a ■ ■ • a f 

Define random variables on Q by 



and denote by the corresponding probability measure on V. : llv' < K~ x • • • x irx x 



lnPjJ(A) lnP^(A^) 



iV AT 
where A = {Ai, A~2, . . . } G £1 

Theorem 2.1 (Asymptotic equipartition property). Let V be a real analytic function grow- 
ing faster than log(l + A 2 ), as |A| — > oo whose equilibrium density p v has q-interval support 
o~ (q > 1). Assume also that p v is generic, that means that p v ^ in the internal points of 
a, p v behaves like square roots near the edges of a, and the function 

v{X) : = 2 J log |/i - \\p v {p)dp - V(X) 

attains its maximum only if A belongs to a. Then for any (3 > the random variables 
converge F^-almost surely to some constant Ep(V). 

Remark 

In the Hermite case, 

^'^)-,nr( 1 + f) + f,( 1 + f)-f-I 

The results of this theorem and Theorem \2.3\ below also hold for the Circular, Laguerre, and 
Jacobi beta ensembles. In particular, 



Ef r = ln(2vr) - lnT ( 1 + £ ) + ^ ( 1 + J 



E%« = ln(2vr) - InT (l + + §V (l + § ) 



^ l ^ fli a 

— — + a In — — at — pmp, 



£^ = ln(27r)-lnr^l + !j + |V ^1 + |J - | - (/i + z/) In 2, 

where ip is the digamma function, 

(io) ^) = ^iogr(z). 

The following is an immediate corollary of Theorem 12.11 

Corollary 2.2. For any (3 > and any e > 0, there exists N Ej p > such that for any 
N > N £ ^, there exists S^(e) C supp(P^) such that for any A G S^(e) we have 

e -( Efi+£ )N < pfiQ) < e -(E p -e)N 

and 

lim / P^{X)dX = 1. 

N ^°° J\ES N (e) 



ENTROPY FOR BETA RANDOM MATRIX ENSEMBLES 5 

2.1.1. The Shannon- McMillan- Breiman Theorem. The interpretation of the constant Ep(V) 
as the entropy of the corresponding random matrix ensemble is by analogy with the Shannon- 
McMillan-Breiman theorem. To illustrate the analogy, we briefly recall the Shannon-McMillan- 
Breiman theorem in the case of the Bernoulli process. 

Let W be the set of sequences of Os and Is, i.e. W = {0,1 } N . Given ro G W let 
ro n = (wi, . . . ,w n ) G {0,1}" be the first n elements of ro, and let S Kn be the cyllinder 
set corresponding to tv n . S Kn is the set of all sequences in W the first n elements of which 
coincide with tv n . Given p G (0, 1), let P p be the Bernoulli measure on W which to a cyllinder 
set assigns the probability 

p(s mn )= P k (i- P r- k , 

where k is the number of ones in ro n . Let 

H(p) = —plnp — (1 — p) ln(l — p) 

be the entropy of the Bernoulli measure. The Shannon-McMillan-Breiman theorem states 
that for any p G (0, 1), we have 

P p (vo G W : lim J^E^A = H ( p \) = L 

Y rnoo n J 

In other words, the theorem states that the random variables _ lnJ p( Sro «) converge to H{p) 
almost surely. 

2.2. Central Limit Theorem. 

Theorem 2.3 (Central Limit Theorem). Under the assumptions of Theorem \2. 11 

Y£:=N-^(lo g pP(\) + NEp(V)) 

converges in distribution as N — ^ oo to the Gaussian random variable with expected value 
and variance 

2 4 v V 2 
where if) is defined in IHU\) . 

Remark 

The results of Theorems \2.1\ and \2.3\ in the special case of Hermite beta ensemble (T3|) were 
first proven by Ionel Popescu in |2"4"] . 
Remark 

Let us write the probability density (QJ] as 

_ e -PH N {\)/2 

where the potential energy term if at is given by 

N 

(11) H N (\) = N^ViXi) - ^ln |Ai - A,|. 

i=l tyj 

Then the Central Limit Theorem result can be reformulated in terms of convergence in dis- 
tribution of 



W N (X) 



N 
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to the Gaussian limit N(0, i — ip' (l + |)), where C^^ is an appropriate centering constant. 
In particular, for the classical beta ensembles, one has 

eft =(\ + if ) - iiv i„ * - i + M-»(i + |) W[ 



ijac 



= !f (AT - 2)AT - iiVlnAf + (-\ In | + i</> ( 1 + § ) ) V. 



cg5 = -^ lllW -(i h f-^(i + f)) W . 



Remark 

It should be noted that both terms in the expression < f77)j /or me potential energy, namely 
N^2 i=1 V(Xi) and Si<i<j<iv^ n — Aj- j , have fluctuations of the order N, (see e.g. [18] , 
[215] ). ili me same time, the sum of these two terms fluctuates on a much smaller order, 
namely N 1 / 2 . The cancellations take place because of the Euler- Lagrange variational equations 
|2Hi|) for the equilibrium measure. As a result, the difference V(A») — 2 I A? ~~ ^ s 

eowa/ (up to a negligible error term) to a constant independent of 1 < i < N. 

It is also important to note that while the fluctuation of the linear statistic ^2f =1 V{\) is 
(asymptotically) Gaussian in the one-cut case, it is non-Gaussian, in general, in the multi- 
cut case (see [23] ). 

For the Hermite beta ensemble, the potential energy H^(X) attains its (unique) minimum 
at the configuration given by the zeroes of the (rescaled) Hermite polynomial of order N, 

namely ^(yfi), where 



[N/2] 



h N (x) = N\J2 



[-l) m {2x 



x N-2m 



m=0 



m\(N - 2m)! 



(see e.g. [201 A.6]). The minimum value of if at is equal to 

(12) N ( N ~ 1) (logiV + 1) - 3 J = - liVlog N-- A N + o(N). 

Note that the inequality + x) > logx implies that 

C% e Z > -N 2 - —NlogN - -N. 

Thus, for a typical configuration the difference between the potential energy and its mini- 
mumal value is proportional to N. 

Similar explicit computations could be done for the other classical beta ensembles. In 
particular, for the Circular beta ensemble, the ground state is unique up to a rotation and is 
given by the vertices of a regular N— gon, and the minimum of the potential energy is equal 
to —^NlogN. For Laguerre and Jacobi beta ensembles, the ground state is unique and is 
given by the zeroes of the corresponding (rescaled) orthogonal polynomial of degree N. Since 
the computations are very similar to the ones in [20} A.6]) we leave the details to the reader. 
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3. Proof of Theorem 12.31 

We first give the proof in the case of the Hermite beta ensemble. Then we quickly extend 
to analytic potentials V in both one-cut and multi-cut cases. 

We start with the Selberg type integral corresponding to the Hermite orthogonal polyno- 
mials [22 17.6.7] 



(13) Z N 03) = (2*) ^— j 11 r(1 + (3/2) . 



which can be rewritten as 



log z»)oj) 4 log 2. - + £ (1 - §)) log ^ - iVlogr(l 4- P72) 

+ ^logr(l + /3j/2). 
J'=l 



Let us represent Y^jLi^°S^{^- + /5j/2) in terms of the Barnes G-function defined in pTJ (see 
also [16] formula (4.181)). It satisfies the functional equation 



G(l + z) = T{z)G(z), G(l) = l, 



and so at the integer points C7 can be represented as G(l + N) = Yif=i^(^ + j)- The 
asymptotic expansion for logG is known (see [T6] formula (4.184)): 



(14) logG(l + z) = ylog2;-^ 2 + ^log27r-^log2 ; + C / (-l) + o(l), 2^00. 



To obtain the representation for Ylf=i 1°S ^(1 + j3j/2) we use the Stirling formula in the form 



(15) \ogT(l + z) = (z+^)\ogz-z+^\og27i + -^- + r(z) \r(z)\ < Cz~ 2 , 
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where the bound for r(z) is uniform for any interval (8, oo) with 8 > 0. Then 

jr log r(i + «/2) = jr ((^ + \) log f - ^ + i log 2^ + -L + ((«) 



- - - - - 6 « 



2 V 2 ; 6J v 2 2 y 6 2 2 V 2 ; 6 12j v 2 (3 J 



= E§ iogr(i + i) + |(i-f) iogr(i + iv) + + f (1 + §)) log 

+ j(l - f ) log2vr - i(| - |) logiV + CW(/3) + o(l) 
J- logG(l +JV) + |(l + f) logr(l + iV) + + f (1 + §)) log I 

(16) + y (1 - f ) log 2tt - 1 (| - |) log JV + C(/3) + o(l). 



Here 



C(/3) = lim CW(/3), 



Af-s>oo 
N 



^C9) = E(^'/2)-fr(;))-^(f-|)7 
i=i 



where 7 denotes the Euler constant and r(z) is defined in (floT) . 
This expansion combined with ( 114)) and (II 5p implies 

v^, n /x Z^ 2 , 3„ „ N , B s , /3N N, B, 

E lo S r ( X + #/ 2 ) = V lGg 2 " 8^ + T (1 + 2 ) ° g 2 " T (1 + 2 ) 

i=i 

+ ^ log 2vr + R(B) log iV + <7(/3) + o(l), 

where 

C(/3) = C(/J) + i(l + |)log2 I r + C'(-l). 
Applying the Selberg formula ( fl3|) . the above representation yields 

log^ er (/3) =f N 2 £ + SI log ^ + Ar(log2vr - lo g r(l + £) - 1(1 + £ 

(17) + R(0)]ogN + C(0) + o(l), Eo = -\. 
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Consider now the random variable 

yHer, P = N -V2 {logp Her,^- x) + N ^ = ^_ Iff^) - \ gZ N (/3) + N E^) . 

The characteristic function of Y^ er,t3 has the form 

$(*) = (Z^ r j dAexp { - §(1 + ^)^(A) - logZ^/J) + tV^f"). 



log$(t) = log ^(0(1 + -£=)) - (1 + log^(^) + (3tVNEf er 



Thus, 

Since according to (|T7j) 

(18) logZ# er (/3) = ^iV 2 £ + ^ logiV + Nf (P) + R{p) \ogN + + o(l) 



with 



/„(/*) = lo g 2vr - logr(l + f ) + f log f - ~(1 + §) 



one can see immediately that 

log^(/3(l + JL)) - (1 + J=) log^(/3) = ty/NtffM - fo(P)) + ^f'(P) + 
v iV v iV ^ 



Hence, using that by definition of i£j rer 

Ef~ = f (p)-l3f (P), 

we get 

$(t) = exp{^f(/3) + o(l)}. 

Thus, we have proved that Y^ er ^ converges in distribution to the Gaussian random variable 
with zero mean and the variance 

m3)=§-^V(i+§)- 

Consider now the case of general one-cut potential V. We write 

n 

(19) H N (X; V) = - JV^y(Ai) + ^log |A, - A, | 

(20) P£(A; 1/) =(Z^C8; v))-i e -^(X;V)/2 | 

(21) Z®{P;V)= [ e-^^' 2 dX. 
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It was proven in [27] (see also [7J) that in the case of analytic one-cut potential, 

(22) log Z N {P- V) =^£[V] + F P {N) + NS[p v ] (| - l) + rp[p v ) + O^' 1 ) 

BN 2 

F,(N)= log Z§-(P)-^S 
S[p v ]=((\ogp v ,p v )-l + \og27r), 

£[V)=-(V,p v ) + (Lp v ,p v ), (Lf)(X) = J\\-fi\f(p)dp 

where p v is the equilibrium density, rp[p] is a smooth function of (3 defined by a contour 
integral with the Stieltjes transform of p v , and (/, g) := J f(x)g(x)dx. Since (122]) is also 
written in the form f fl8|) with £ replaced by £[V], and /o, C((3) replaced by 

fviP) = hiP) + S[p v ] (| - l), C[V) = C(J3) + rp[p v ], 
by the same argument we conclude that the CLT is valid for 

Y$ = ATV^iogP^A; V) + NE P [V\) = N^ 2 ^- ^H N (X; V) - \ogZ N (P; V) + NEp[vfj, 
where 

Ef,\y\ = Ef« - S\p v ], 
and the variance is the same as in the Gaussian case. 

We are left to consider the multi-cut case. Since the analog of (1221) obtained in [28] in the 
multi-cut case is more complicated, we need first some extra definitions. 

Denote 

i 



(23) a = (J cr Q , cr a = [a a ,b a ], p a = / p^(X)dX, pi := l CJa p i 

„ 1 J Gci 



p a yA)UA, p a := ± a<x f 

a=l J<T * 

Define the operators £, C and £ 

(24) Cf = l a L[fl a ], C a f := la a L[fl«J, £ ■= < =1 L, £:=£-£, 

on the set of the functions 

H = ®l =x L x [a a \. 

Note that the topology of H is not important below. 
For each interval a a we also define the operator 



(25) D a = \D a + DH D a h(X)- 1 f h\»W{p-a a )(b a -p)dp 



2 V an ^ Ja a (X-fiW(X-a a )(b a -X) 

and set 

(26) D := ® q a=l D a . 
Let Q be a positive definite q x q matrix of the form 

(27) Q = {Q aa >} q a , a , =1 , Qaa' = -(^ (a) ,^ (Q,) ), 

where the function ^ a '(X) is a unique solution of the system of equations 

(28) - (£ip {a) ) a > = <W, a' = l,...,g. 
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0W(A) = p Q (A)X a - 1 / 2 (A)l (T , xy 2 (X) =$*(j[(z- a a ){z - b a ) 
where p a is a polynomial of degree q — 1. Set 

(29) I[h] = (h[h},..., I q [h]), I a [h) :=£Q^(/^ K) ). 
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1/2 



z=\+iO 



and define a quasi-periodic in N function 0(/3; {iVp}): 





(30) 003; {iVp}) := £ exp { - | (V An, An) + (| - l) (An, Ipog? v ])}, 

niH \-n q =no 

g 

{Nfl} = ({N t x 1 },...,{NfM q }), (An) a = n a -{Nfi a }, n = ^{iVpJ, 

c*=l 

with Q of (J27D, {pJLi of ([23]), J[/i] of (J29J), and logp = (logpi, . . . , logp,). 

Using the above definitions we can write now the analog of (122]) obtained in [28] : 

log 4 6) [/3, V] =-£*\y\ + F P {N) + iV(| - + (q - 1) (J2(0) log iV + C(/3)) 



(31) 



9 1 — ~ 

+ I>^aV] + -ft(yS) log/io) - -logdet(l - DC) 



a=l 



+ |(f - l) 2 (^(l " DC)~ 1 V,V} + log0(/3; {iVu}) + 0(iV" K ), « > 0, 



where p a are defined in ( 123]) . r j g[/i Q , 1 p Q ,] for each a a is the same as in ( [22]) . Fp(n), R((3), 
and C(/3) are defined in ([17]) . and det means the Fredholm determinant of _D£ on cr. The 
non positive measures (z/i, . . . , u q ) =: 17 have the form 



( u h) :=-(h(b a ) + h(a a )) - — I h^d X — 

V ' 4 l 1 a) K a,> 2tt I y/(X-a Q )(b a -X) 2 



+ -(D a log P a ,h). 



with 



2vrpr(A) 



a/(A - a a )(b a - A) 

It is easy to see that in the multi-cut case the expression (]3T]) still can be written in the form 
( 11%]) with £ replaced by £[V], -R(/3) replaced by qR((3), and / , C(/3) replaced by 

M/3) =/o03) + % v ](f-i). 

9 1 _~ 

CJF] =gC(/3) + ^(^[^VJ + R(P)\o gf i a ) - -logdet(l - 



+ |(f - l) 2 (^(l " DC)-%V) + log0(/3; {iVp}). 



Since the structure of 0(/3; {iVp}) (see ([30]) ) guarantees that 



lo g e(/3;{jV//}) -lo g e(/3(l + 



< 



N 
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by the same argument we conclude that the CLT is valid for 

Y$ = iV- 1 /2(logPjJ(A; V) + NE P [V}) = N~ l ' 2 (^H N {\; V) - \ogZ N {fr V) + NE P [V\), 
where 

E P \V] = E^-S[p v ], 

and the variance is the same as in the Gaussian case. 
Theorem is proven. 
Remark 

As we have mentioned above, the result of Theorem \2.'J\ also holds for Circular, Laguerre, 
and Jacobi beta ensembles. The proofs are very similar to the one in the Hermite case. The 
arguments are based on the explicit formulas for the partition functions given by the Selberg 
integrals [20J. 



Z N {?) =(2rr) ni + f}/2)N , 

TT r(l + U + 1 W2)r(/i + jP/2)T(v + j/3/2) 

N {P) f} r(i + /3/2)r(/i + ^ + (iv + i-i)/3/2) ' 

The details are left to the reader. 
To prove Theorem 12.11 it suffices to observe that by the Chebyshev inequality 

Prob{\X p N + E p \ >e} = Prob{\Y^\ > VNe} < ($(1) + $(-l)) e -^ e < Ce~^ £ 
and apply the Borel-Cantelli lemma. 
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